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Abstract 

Background: Differential cross sections for the (p, n) and (n,p) reactions on 90 Zr over the interval 
of — 50 MeV excitation energy were used to determine the corresponding GT strengths, and 
the resulting quenching factor ~ 0.9 with respect to the Ikeda sum rule. In this procedure the 
contribution of the isovector spin monopole (IVSM) strength was subtracted from the total strength 
without taking into account the interference between the GT and the IVSM modes. 

Purpose: To determine the quantitative effect of the IVSM excitation mode on the L = strength 
in charge-exchange reactions on several closed-shell nuclei and the Sn isotopic chain. 

Method: The fully consistent relativistic Hartrcc-Bogoliubov (RHB) model + proton-neutron 
quasiparticle random phase approximation (pn-RQRPA) is employed in the calculation of transition 
strength in the f3~ and /3 + channels. 

Results: The inclusion of the higher-order terms, that include the effect of finite momentum 
transfer, in the transition operator shifts a portion of the strength to the high-energy region above 
the GT resonance. The total strength is slightly enhanced in nuclei with small neutron-to-proton 
ratio but remains unchanged with increasing neutron excess. 

Conclusions: Terms that include momentum transfer in the transition operator act mostly to 
shift the strength to high excitation energies, but hardly affect the total strength. Based on the 
strength obtained using the full L = transition operator in the pn-RQRPA calculation, we have 
estimated the impact of the IVSM on the strength measured in the charge-exchagne reactions on 
90 Zr and found that the data are consistent with the Ikeda sum rule. 

PACS numbers: 21.60.Jz, 24.30.Cz, 25.40.Kv 



I. INTRODUCTION 



Collective spin and isospin excitations in atomic nuclei have been the subject of many 
experimental and theoretical studies (for an extensive review see Ref. |l|). Of particular 
interest is the Gamow- Teller (GT) resonance, a collective oscillation of neutrons that coher- 
ently change the direction of their spins and isospins without changing their orbital motion. 
A detailed knowledge of the GT strength distribution is essential for the understanding of 
nuclear beta-decay and weak processes in stars [4J]. It was first predicted in 1963 [2| and 
observed in (p,n) reactions a decade later [3j. Further measurements identified a problem 
that is still actively discussed, and that is the quenching of the Gamow- Teller strength. In 
numerous experiments across the whole nuclear chart only around 60% of the strength pre- 
dicted by the model- independent Ikeda sum rule had been observed |l| . These experiments, 
however, were only able to measure the strength up to the excitation energy of the giant 
resonance. A consistent analysis of (p, n) and (n,p) reaction data from 90 Zr over a much 
wider range of excitation energies, concluded that the GT strength is actually quenched 
by approximately 10% [5|, |6|. Theoretical models systematically overestimate the transition 
strength compared to the measured values |7|. This effect was attributed to two possible 
processes: (i) coupling of the GT mode to A-isobar nucleon-hole (A — h) configurations, 
and (ii) second-order configuration mixing. It has been shown that the former process is 
responsible for only a small fraction of the quenching [8|, l9j|, leaving the latter as the major 
mechanism for shifting the GT strength to higher energies 10] . 



The spin-isospin operator structure of the (p, n) probe is similar to that of the Gamow- 
Teller (GT) operator [8j. However, they become comparable only if the GT cross section is 
measured at very small momentum transfer q. In the (p, n) reaction this condition can be 
met only for zero degree scattering, small excitation energies, and high bombarding energies. 
Extraction of the L = strength at high excitation energies, where higher-multipole response 



dominates, is very difficult ll|, |l2|. Nevertheless, recent experiments reported data on the 



Gamow- Teller response in 90 Zr, in both (3~ 5] and (3 + channels [6j, up to 50 MeV excitation 
energies. Therefore, in the total strength, contributions from higher-order terms in the 
expansion in q appear, the first of which is the isovector spin monopole (IVSM) mode. This 
mode, with the transition operator r 2 crr, represents a collective excitation of the nucleus with 
quantum numbers J 7T = 1 + ,L = 0,S = 1 and T = 1. Even though the first observation 



of the isovector spin monopole mode was reported in 1983 13J, a quantitative analysis 
and determination of the strength remain difficult 14H16J]. The unknown IVSM strength 
also introduces uncertainties in the measurement of the Gamow- Teller strength. To obtain 
precise quantitative data on the total GT strength, and indirectly on the quenching, the 
contribution of the isovector spin monopole mode and of higher-order terms in momentum 
transfer, must be subtracted from the measured strength. Although interference occurs 
between the GT and IVSM modes, these contributions are usually subtracted incoherently 
from the spectrum because the distribution of GT strength in the IVSM resonance region 



is unknown [6|. Related to new measurements of nuclear response in unstable nuclei 17] . 
the correct treatment of effects that influence the extraction of the Gamow- Teller strength 
is all the more important. 

Not much theoretical work has been reported on the IVSM strength so far. Isovector spin 
excitations with angular momentum L = 0,1, and 2 have been studied in Ref. 18| employing 
the Skyrme SIII Hartree-Fock model and the random-phase approximation (RPA) with a 
schematic residual p-h interaction. The Skyrme functionals SGII and SIII were used in a 
self-consistent HF + Tamm-Dancoff approximation (TDA) study that focused only on the 
IVSM mode in 48 Ca, 90 Zr and 208 Pb [19J. In both cases the IVSM mode was identified at 
excitation energies between 20 MeV and 60 MeV with respect to the ground state energy 
of the parent nucleus. A non-energy-weighted sum rule was devised that involves particle 
numbers and radii, and the calculated strengths were shown to be consistent with values 
obtained from ground-state densities. The effect of the isovector spin monopole mode and 
higher-order terms in momentum transfer on the quenching of GT strength has not been 
investigated so far. 

In this study we explore the IVSM mode of excitation within a self-consistent micro- 
scopic theory, and analyze the effect of momentum transfer on spin-isospin excitations. The 
relativistic Hartree-Bogoliubov (RHB) + proton-neutron relativistic quasiparticle random 
phase approximation (pn-RQRPA) framework is employed for the calculation of the IVSM 
strength. This framework, based on the covariant energy density functional theory, is a 
charge-exchange extension of the relativistic quasiparticle RPA formulated in the canonical 
basis of the RHB model [20J . The RHB + pn-RQRPA have already been successfully app lied 
to the analysis of the Fermi and the Gamow- Teller response 21[ , /3-decay half-lives 22|, |23j , 
neutrino- nucleus cross-sections [2J|, total muon capture rates 25], and electron capture 



rates [261 ] . 



The IVSM strength is calculated for the closed-shell nuclei 48 Ca, 90 Zr and 208 Pb. The 
isovector spin monopole operator is also considered in the context of an expansion of the 
transition operator with respect to the momentum transfer in the reaction. The impact of 
the momentum transfer on the total L = strength and its distribution is examined in the 
Sn isotopic chain between A = 100 and A = 150. Sec. [Ill introduces the formalism, and Sec. 
Ill II presents the results and discussion. Sec. IIVI contains a short summary and concluding 
remarks. 



II. THEORETICAL FORMALISM 

The relativistic quasiparticle random phase approximation (RQRPA) was formulated 
in the canonical single-nucleon basis of the RHB model in Ref. [20] and extended to the 
description of charge-exchange excitations (pn- RQRPA) in Ref. [21] . The RHB + RQRPA 
model is fully self-consistent: in the particle-hole channel, effective Lagrangians with density- 
dependent meson-nucleon couplings are employed, a nd p airing correlations are described by 



the pairing part of the finite range Gogny interaction [27] . In both the ph and pp channels the 
same interactions are used in the RHB equations that determine the canonical quasiparticle 
basis, and in the matrix equations of the RQRPA. This is very important because the energy 
weighted sum rules are fulfilled only if the pairing interaction is consistently included both 
in the static RHB and in the dynamical RQRPA calculation. In the present work all cal- 
culations are performed using one of the most accurate meson-exchange density- dependent 



relativistic mean-field interactions in the ph channel: DD-ME2 28]. 

Transitions between the + ground state of a spherical even-even parent nucleus and the 
J n excited state of the corresponding odd-odd daughter nucleus are induced by a charge- 
exchange operator T JM . Taking into account the rotational invariance of the nuclear system, 
the quasiparticle pairs are coupled to good angular momentum and the matrix equations of 
the pn-RQRPA read: 



Ml... (i) 





The matrices A and B are defined in the canonical basis 29] 



J ^pn,p'n' rL ppi°nn' \ £1 nn'°PP l 



+ (u p V n U p/ V n ' + V p U n V p iU n <) V pn , np , + (u p U n U p >U n > + V p V n Vp>V n >) V^ p , n , 
B Lp'n' = (-l) j r'- 3n ' +J {u p V n V p m n , + V p U n U p ,V n ,) V^n' 



-pn,p'n' \ ~J \"-p~n»p ~tj • -p»Ti«ji - n j • pp < 

yppj 
pnp' 



(u p u n v p iv n i + v p v n u p <u n i) V?L, n , . (2) 



Here p, p', and n, n' denote proton and neutron quasiparticle canonical states, respectively, 
V ph is the proton-neutron particle-hole residual interaction, and V pp is the corresponding 
particle-particle interaction. The canonical basis diagonalizes the density matrix, and the 
occupation amplitudes v P>n are the corresponding eigenvalues. However, the canonical basis 
does not diagonalize the Dirac single-nucleon mean-field Hamiltonian ho and the pairing 
field A, and therefore the off-diagonal matrix elements H^, and H ppl appear in Eq. fl2]): 

Hll, = (u K u K > - v K v Kl )h KK , - (u K v K > + v K u Kl )/\ KK i , (3) 

For each energy E\, X XJ and Y XJ in Eq. ([T]) denote the corresponding forward- and 
backward-going QRPA amplitudes, respectively. The total strength for the transition be- 
tween the ground state of the even-even (N,Z) nucleus and the excited state of the odd-odd 
(N+1,Z-1) or (N-1,Z+1) nucleus, induced by the operator T JM , reads 



Bfj 



J2 < P\\T J \\n > {X^u p v n + (-l) J Y p x n J v p u n ) 



pn 



(4) 



The discrete strength distribution is folded by the Lorentzian function 

r 

71 



i?(E) ± = y b±- 2 — 2 , . (5) 



In the present calculation the width of the Lorentzian function is T = 1 MeV. 

The spin-isospin interaction terms are generated by p-and 7r-meson exchange. Because of 
parity conservation, the one-pion direct contribution vanishes in the mean-field calculation 
of a nuclear ground state. Its inclusion is important, however, in calculations of excitations 
that involve spin and isospin degrees of freedom. The particle-hole residual interaction in 
the pn-RQRPA is derived from the Lagrangian density 

4 n t = -9 P ^P^ - h-^rfQfi^ . ( 6 ) 

TTljr- 



Vectors in isospin space are denoted by arrows, and boldface symbols indicate vectors in 
ordinary three-dimensional space. 

The coupling between the p-meson and the nucleon is assumed to be a function of the 



vector density p v = y/J^J^, with j^ = ^j^. In Ref. [30| it has been shown that the explicit 
density dependence of the meson-nucleon couplings introduces additional rearrangement 
terms in the residual two-body interaction of the RRPA, and that their contribution is 
essential for a quantitative description of excited states. However, since the rearrangement 
terms include the corresponding isoscalar ground-state densities, it is easy to see that they 
are absent in the charge exchange channel. For the p-meson coupling the functional form 
used in the DD-ME2 density-dependent effective interaction [28( reads 

g P {Pv) = g P {Psat)e- a ^-v , (7) 

where x = p v /p S at, and p sa t denotes the saturation vector density in symmetric nuclear 
matter. For the pseudovector pion-nucleon coupling the standard parameters are used (see 

Ref. [311), 

f 2 
m„ = 138.0 MeV -^- = 0.08 . (8) 

The derivative type of the pion-nucleon coupling necessitates the inclusion of a zero-range 
Landau-Migdal term, which accounts for the contact part of the nucleon-nucleon interaction 

V S7 r=g'(—) fi^Si • S 2 (J(n - r 2 ) , (9) 



where 




(10) 



and the parameter g' DD _ ME2 = 0.52 is adjusted to reproduce the GTR excitation energy in 
208 Pb 

The pn-RQRPA model is fully consistent: the same interactions, both in the particle- hole 
and particle-particle channels, are used in the RHB equation that determines the canonical 
quasiparticle basis, and in the pn-RQRPA Eq. ([[]). In both channels the same strength 
parameters of the interactions are used in the RHB and RQRPA calculations. With respect 
to the RHB calculation of the ground state of an even-even nucleus, the charge-exchange 
channel includes the additional one-pion exchange contribution. 



The two-quasiparticle configuration space includes states with both nucleons in the dis- 
crete bound levels, states with one bound nucleon and one nucleon in the continuum, and 
also states with both nucleons in the continuum. In addition to configurations built from 
two-quasiparticle states of positive energy, the RQRPA configuration space contains pair- 
configurations formed from fully or partially occupied states of positive energy and empty 
negative-energy states from the Dirac sea. The inclusion of configurations built from occu- 
pied positive-energy states and empty negative-energy states is essential for the consistency 
of the model 21] . 



In the pp-channel of the RHB model a phenomenological pairing interaction is used, the 
pairing part of the Gogny force, 

V pp (l,2) = J2 e- {{ri - r2)/ ^ )2 {Wi + BiP a -H Z P T -MiP a P T ), (11) 

i=l,2 



with the set D1S 32| for the parameters fa, Wi, B>i, Hi and M, (i = 1,2). This force has been 
very carefully adjusted to the pairing properties of finite nuclei all over the periodic table. 
In particular, the basic advantage of the Gogny force is the finite range, which automatically 
guarantees a proper cut-off in the momentum space. The same Gogny interaction is also 
used in the T = 1 pairing channel of the pn-RQRPA. 

III. RESULTS 

In the first part the calculated strength distributions for the isovector spin monopole 
transition operator in 90 Zr and 10 °- 150 Sn are analyzed. We then calculate the L = strength 
and show the effect of momentum transfer on the strength distribution and the total strength 
in the tin isotopic chain. 

A. Isovector spin monopole strength 

The isovector spin monopole (IVSM) operator reads 

T /VSA/ = ^ r 2 Sr± _ (12) 

i=l 

In Fig. [T]we display the IVSM strength for 90 Zr up to 70 MeV excitation energy. Because 
of the structure of the IVSM operator of Eq. (IT2"|) . the dominant feature of the spectrum is 

8 



8000 



> 6000 - 90 



<*" 4000 
E 

H — m 

cr 2000 - 



30 40 

E [MeV] 




i 

--- with Ofi co 
— without Oh co 



60 



70 



FIG. 1. (color online) The pn-RQRPA strength distribution of isovector spin monopole states in 
90 Zr. The dashed curve denotes the total strength with 0/iu; transitions included, whereas the solid 
curve corresponds to the strength for which only 2hu) and higher excitations are included in the 
configuration space, for the j3~ and the f3 + channels in the upper and lower panels, respectively. 

a strong peak at the position of the Gamow- Teller (GT) resonance. However, unlike the GT 
operator that excites only Ohu transitions, the isovector spin monopole operator can also 
excite 2hu transitions. These relatively weak transitions at excitation energies above the 
giant resonance contribute a significant portion of the total strength, because even though 
each individual transition is weak, their number is large. The 2Hcu transitions are also 
responsible for all the strength above the GT resonance in the /3 + channel. The amount of 
strength above the resonance in the two channels is comparable, so that the value of the 
sum rule is mostly determined by the strength contained in the resonances. 



The non-energy-weighted sum rule for the IVSM transition strength reads 18] 



O— — O-l 



(13) 



3[N(r^ n -Z(r^ p 
Using the values of (r 4 ) nv that correspond to the RHB self-consistent ground-state solution, 

9 



<n,p 



we have verified that the sum rule Eq. f)13p is satisfied in our calculation (cf. Table H]). 

TABLE I. Integrated strengths of the isovector spin monopole operator in 48 Ca, 90 Zr and 208 Pb 
nuclei. Proton and neutron radii correspond to the RHB self-consistent ground solution. All values 
are given in units of fm 4 . 
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Since the IVSM strength distribution is dominated by the Ohu components, to study 
the behavior of the IVSM mode it is convenient to exclude these transitions. This can be 
done by including only 2hu and higher configurations in the QRPA basis. The solid curve 
shown in Fig. [TJ corresponds to the strength for which only 2hw and higher excitations 
are included in the configuration space. In Table [Til we compare the energy centroids of 

TABLE II. Energy centroids of the isovector spin monopole strength in the f3~ and the f3 + channel, 
for the nuclei 48 Ca, 90 Zr and 208 Pb. The values calculated with the relativistic functional DD-ME2 
are compared to those obtained with the Skyrme functionals SGII and SIII 



19] 
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the IVSM distributions in the j3~ and the /3 + channels (excluding Ofkv configurations) , 
with the corresponding values obtained with two Skyrme functionals: SGII and SIII |19| . 
The agreement is very good in the (3~ channel, whereas in the j3 + channel the centroids 
calculated with the relativistic functional DD-ME2 are found to be few MeV higher than 
those predicted by the two Skyrme functionals. Both calculations predict a decreasing of 
the energy centroids with increasing mass, in agreement with the results of Ref. 18 1. 



10 



50 



40 



^30 
> 

CD 



LU 



20 



10 







¥ , I , I , I , I I I I I I I I I I I , 

■ 


1 ' 1 ' 


1 ' 


■ 
■ 
■ 


• GT 
■ IVSM 


- 


■ 

■ ■ 


■ 

■ ■ 

■ 


■ ■ 


" Sn 

1,1,1,1,1,1,1,1,1,1,1,1,1, 



100 104 108 112 116 120 124 128 132 136 140 144 148 

A 



FIG. 2. (color online) Energy centroids of the isovector spin monopole strength in tin isotopes 
(excluding the Ohuj transitions), in comparison with the corresponding Gamow- Teller centroids. 

For the Sn isotopic chain with neutron-to-proton ratio ranging from N/Z = 1 to N/Z = 2, 
in Fig. [2] we compare the centroids of the Gamow- Teller and the IVSM strength, considering 
only the 2fiw and higher transitions for the latter. It appears that the IVSM strength function 
has a somewhat more pronounced mass and/or isospin dependence. In the lighter Sn isotopes 
the IVSM centroids are found just below 50 MeV, or approximately 25 MeV above the GT 
centroids. At mass 150 the IVSM energy centroid rapidly approaches 20 MeV, only 13 MeV 
above the corresponding GT centroid. 

B. The L=0 strength 

The L = strength obtained in charge-exchange reactions corresponds to the squared 
matrix element of the L = operator 



T(±) = JL=o(qr)^r± 



(14) 



11 



where q is the momentum transfer. In the long wavelength limit (i.e. q — > 0) the spherical 
Bessel function can be approximated by 

j ( qr )Kl-'L- + ... , (15) 

and usually only the first term is retained. However, if the momentum transfer is not 
negligible then, together with the Gamow- Teller operator, the isovector spin monopole term 
has to be taken into account 



Q 



2 



0(±) = Sr ± - jr 2 12r ± . (16) 

For a (p, n) reaction the total energy and momentum of a proton with kinetic energy T 



E P = T + m p , p p = ^jEl - m 2 , (17) 

and for the outgoing neutron: 

E n = E p - E x , Pn = sJWl - m 2 , (18) 

where E x is the excitation energy of the nucleus with respect to the ground state of the 
target (parent) nucleus. The momentum transfer is defined as 



\q\ = \Pp-Pn\ = y P% + Pl- 2 PpPn COS $, (19) 

where d denotes the angle between the momenta of the incoming and outgoing particles. 
Assuming the cross section is measured at forward angles, one can set •§ ~ 0°, and obtain a 
simple expression for the momentum transfer: 

q = \P P ~Pn\- (20) 

Using Eqs. (IT7j) - fl20l) . one notices that the momentum transfer depends linearly on the 
excitation energy of the nucleus, and has a 1/(1 + T/m) 2 dependence on the kinetic energy 
of the incoming proton, as shown in Fig. [31 This makes the effect of higher-order terms in 
the expansion Eq. (1151) more pronounced for higher excitation energies and lower incoming 
energies. 

Fig. H] shows a comparison between the L = strengths in 90 Zr, calculated with the 
full operator Eq. (1141) . the Gamow- Teller operator, and the g 2 -order operator Eq. (1161) . In 
the upper panels we display the strength distributions in the (3~ channel. At excitation 
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FIG. 3. Momentum transfer for a (p,n) reaction, calculated using Eqs. (|17|) - (|2U|) . In the upper 
panel q is plotted for a constant kinetic energy of the incoming proton: T = 300 MeV, and in the 
lower panel for constant excitation energy of the nucleus of E x = 50 MeV. 

energies below 30 MeV, shown in the left panel, momentum transfer is rather small. The 
only significant contribution to the total strength below the GT resonance comes from the 
Gamow- Teller term, while in the region of the resonance the IVSM term in the operator 
actually reduces the strength calculated with the Gamow- Teller operator. 

The largest contribution to the strength of the Ohu part of the GT + IVSM operator 
comes from the orbits around the Fermi surface. In the harmonic oscillator basis the mean 
value of the r 2 is equal for all orbits in a major shell, hence the following proportionality 
relation is obtained: 



O(0M = J2 



i -j 



1 



q 2 r 2 
~6~ 



Sr ± 



j > a { dj 



(l - ^ (r 2 )) 53 (i |Sr±l i> a, + a i( (21 



where the (r 2 ) denotes the mean value of the r 2 operator in the major shell. This relation 
implies that the value of the GT + IVSM matrix element will always be lower than in the 
case of the Gamow- Teller operator. The reduction will be greater with increasing momentum 
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FIG. 4. (color online) Comparison of the pn-RQRPA strengths obtained with the Gamow- Teller 
operator (dashed), the GT + IVSM operator Eq. (fl"6j) (full), and the full L = operator Eq. (fl~4"|) 
(dash-dotted) in 90 Zr. The upper panels display the strength in the j3~ channel, whereas the 
strength in the /3 + channel are shown in the lower panels. Different scales are used for the region 
of low excitation energy below 30 MeV (left), and excitation energies in the interval 20 to 70 MeV 
(right). 



transfer, i.e. with increasing excitation energy with respect to the ground state of the parent 
nucleus (see Fig. HJ). 

As already shown in Fig. [TJ the strength at high excitation energies originates from 
2fijjj transitions. For instance, the peak at 32.5 MeV predominantly corresponds to the 
^1/5/2 - > ^2/7/2 transition. The matrix elements of the GT operator are small in this case, 
and the IVSM term of the operator dominates. In the expansion of the spherical Bessel 
function jo(qr) successive terms have alternating signs, so that the next term reduces the 
strength of the isovector spin monopole mode. This is particularly visible above 40 MeV 
excitation energy as the next term is proportional to g 4 /120. 

14 



The corresponding strength in the f3 + channel is plotted in the lower panel of Fig. HI In the 
low-energy region the behavior is similar to that of the (3~ channel, the contribution of terms 
with finite momentum transfer slightly reduces the strength of the GT resonance. However, 
this reduction is smaller due to the lower momentum transfer involved, corresponding to 
the lower excitation energy. At excitation energies above 20 MeV, the strength is strongly 
suppressed by higher-order terms in expansion in Eq. (1151) . in contrast with the (3~ channel. 
Thus instead of shifting the strength to higher energies, in the /3 + channel finite momentum 
transfer simply reduces the total strength. The enhancement and reduction of the strength 
in the (3~ and (3 + channels, respectively, was observed in the nuclei studied in this work. 

The shift of the strength to higher energies is further analyzed for the Sn isotopic chain 
and illustrated in Fig. |5j where we show the effect of finite momentum transfer on the 
strength distribution and on the total strength in the (3~ channel. The energy E$ 5 % below 
which one finds 95% of the calculated Gamow- Teller strength, is defined by the relation 

" ^ E,* (G £ ( (22) 

E 9 z,% ranges from > 30 MeV in the lightest tin isotopes, to ~ 13 MeV in 150 Sn. For each 
isotope we calculate the ratio of the total strength (including the effect of finite momentum 
transfer) below E 95 % and the total strength, 

v ~ EMX) ' (23) 

where X denotes the operators O and T introduced in Eqs. flT6|) and (TH1) . respectively. 

Values of r\ < 0.95 indicate that the higher-order terms in the operator shift a portion of the 

strength from the resonance to higher excitation energies. To see the effect of momentum 

transfer on the total strength, we also plot the ratio of the total strength calculated with 

the operators O Eq. (fl6|) and T Eq. (TH1) . and the total Gamow- Teller strength 

= luBjjX) 
*> Y,iBi{GTY l } 



For Sn isotopes with masses in the interval 100 < A < 150 the ratios defined in Eqs. (123]) 
and (HM|) . are plotted in Fig. [5j Open symbols denote results obtained with the GT+IVSM 
operator Eq. ffT6l) . One notices that for light isotopes a rather large amount of strength is 
found above the Gamow- Teller resonance. The fraction of strength found at lower energies 
increases with the addition of neutrons to 77 = 0.8 for 116 Sn. From A = 116 the ratio rj is a 
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FIG. 5. (color online) Ratios r\ and (, denned in Eqs. ([23]) and ([23]) . respectively, for the Sn isotopes. 
Open symbols denote ratios calculated with the GT+IVSM operator Eq. ([16]) . filled symbols are 
for those obtained using the full L = operator Eq. (TH]) . 

linear function of the mass, and reaches the value of rj = 0.9 for 150 Sn. The ratio of the total 
strengths ( is considerably larger than 1 for lighter isotopes, but rapidly converges to £ = 1 
with the addition of neutrons. With filled symbols we denote results obtained using the full 
L = operator Eq. (fT4T) . The prominent feature is that, with respect to the GT+IVSM 
operator, the ratio rj increases and ( decreases for all isotopes. Since the next term in the 
expansion Eq. ([15]) is proportional to q 4 , a relatively large momentum transfer is necessary 
for an effect to be noticeable (see also the right panels in Figs. S]andE])- Therefore, this term 
in the expansion does not affect the resonance but reduces the strength at high energies, 
and in this way reduces the total strength and increases the fraction of the strength below 
the GT resonance. Because the energy centroid of the IVSM is higher in isotopes with a 
low number of excess neutrons, with correspondingly large momentum transfer, higher-order 
terms have a more pronounced effect in Sn isotopes with A rj 100. The largest differences of 
the ratios rj and ( with respect to those calculated with the GT+IVSM operator Eq. ([TrJ]) . 
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144 Sn (lower panels), calculated with the Gamow- Teller operator, the GT+IVSM operator Eq. (|16p 
and the full L = operator Eq. (|14p . Note the different scales used for the resonance region and 
the region of high-excitation energies. 



are found in the lightest Sn nuclei. 

This effect is further illustrated in Fig. [6] where we compare the strength calculated using 
the Gamow- Teller operator, with those obtained using the GT+IVSM operator Eq. (|T6|) 
and the full L = operator Eq. ffT4l . in two Sn isotopes: 100 Sn and 144 Sn. In the former 
the resonance is at 23 MeV excitation energy with respect to the ground state of the parent 
nucleus, and the corresponding momentum transfer is q = 0.181 fm _1 for the incoming 
proton kinetic energy T = 300 MeV. The inclusion of the IVSM term reduces the strength 
of the resonance by approximately 30%. In 144 Sn the resonance is at 11 MeV, with the 
corresponding momentum transfer q = 0.086 fm^ 1 . Since the square of momentum transfer 
appears in the IVSM operator, in this case the effect on the resonance is significantly smaller. 
It is important to note that, even though the relative reduction of the resonance is more 
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FIG. 7. (color online) f3~ pn-RQRPA strength located at energies above £"95% in Sn isotopes with 
100 < A < 150. The GT operator (circles), the GT + IVSM operator defined in Eq. (|16p (squares), 
and the full L = operator Eq. (|14p (diamonds), have been used in the calculation of the strength 
distributions. 

pronounced in the lighter isotope, more strength is actually subtracted from the resonance 
in 144 Sn: B{GTJ) - B{TJ) = 8.0, than in 100 Sn: B{GTJ) - B(T_) = 3.7. In Fig. we plot 
the strength that is calculated at energies above £"95% for the Sn isotopic chain. One notices 
that the difference between the strength obtained using the operators defined in Eqs. (fT4"j) 
or fllpp . and the GT strength is practically constant. Because of this in light isotopes there 
is enough additional strength to overcome the reduction of the resonance and even increase 
the total strength. In heavy isotopes the strength at high energies mostly compensates for 
the strength lost in the resonance but does not increase the total strength. 

In Table UTTl we display the values of the ratio ( defined in Eq. (l24"j) . in both the (3~ and the 
(3 + channel for three representative nuclei: 48 Ca, 90 Zr and 208 Pb. The momentum transfer 
corresponds to a kinetic energy of 300 MeV for the incoming proton (neutron). The results 
for 90 Zr are particularly important because of the recent analysis of both (p,n) and (n,p) 
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data 



BE 



6], that determined the GT quenching factor 



gGT _ gGT 

Q = f {N _ Y) = °- 88 ± °- 06 • ^ 



In the (3~ channel the L = strength was measured up to 50 MeV excitation energy: Si- = 
33.5 ±0.6(stat.) ±0.4(MD) ±4.7(<3"gt)- Employing a distorted wave impulse approximation 
(DWIA) model to estimate the contribution of the isovector spin monopole strength, the 
assumption was made that the complete IVSM strength is concentrated in a single state at 
35 MeV. The estimated contribution was then subtracted from the measured strength, and 
the value of the total Gamow- Teller strength was determined: S9- = 29.3 ± 0.5(stat.) ± 
0.4(MD) ± 0.9(IVSM) ± 4.7(<x GT ). This means that the IVSM contribution enhances the 
total strength by approximately 15%. Our results using the full L = operator show (cf. 
Table IIIIj) . that the total strength in the /3~ channel is not modified by the inclusion of 
higher order terms, and their only effect is to shift part of the strength to energies above 
the resonance. Therefore, the result of the present calculation implies S9I = Sir . 

The experimental value of the L = strength in the (3 + channel was determined in Ref . [6j : 
S^=° = 5.4±0.4(stat.)±0.3(MD)±0.9(a G r). After subtracting the IVSM strength, the value 
of the GT+ strength was deduced: S$? = 2.9±0.4(stat.)±0.3(MD)±0.3(IVSM)±0.5(<j GT ). 
The present calculation indicates, however, that the total GT strength in the /3 + channel is 
actually reduced by ~ 15% by the inclusion of finite momentum-transfer terms. Therefore, 
using the total measured L = strength, we deduce the GT strength in the j3 + channel: 
S9+ = 6.3. The deduced value for the Ikeda sum rule: S9- — S9+ = 33.5 — 6.3 = 27.2, is 
consistent with the quenching factor extracted from data in Ref. 6j. However, in the (3~ 
channel data were only obtained below 50 MeV excitation energy, whereas our calculation 
predicts that approximately 6% of the total strength is located above this energy. Assuming 
that the measured strength actually corresponds to only 94% of the total strength, we obtain 
S9- = 35.6 and, therefore, the value of the sum rule: S9- — S9+ = 29.3. Considering the 
experimental uncertainty, in particular the one originating from the Gamow- Teller unit cross 
section, this result may indicate that no quenching of the experimental strength with respect 
to the Ikeda sum rule occurs. We note that arguments for this conclusion were already put 



forward from the point of view of the shell model 33] 
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TABLE III. The ratio ( defined in Eq. ([MD, for the j3~ and /?+ channels in 48 Ca, 90 Zr and 208 Pb. 
The second and third columns display results calculated with the GT+IVSM operator Eq. ([To]) . 
In the last two columns we show the results obtained with the full L = operator Eq. (|14j) . The 
momentum transfer corresponds to a kinetic energy of 300 MeV for the incoming nucleon. 
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IV. 



CONCLUSION AND OUTLOOK 



An accurate determination of Gamow- Teller strength remains a challenge for charge- 
exchange reaction experiments. With the progress of experiments that can provide data 
on the nuclear response at high excitation energies, the effect of finite momentum transfer 
must be taken into account. In this work the L = strength has been analyzed in the Sn 
isotopic chain, 48 Ca, 90 Zr and 208 Pb. Employing the RHB + pn-RQRPA framework, we 
have compared strength functions calculated using the GT operator, the GT plus isovector 
spin monopole mode term, and the operator that contains the full momentum transfer 
dependence. 

The transition strength for the pure isovector spin monopole operator has been calculated 
for 48 Ca, 90 Zr and 208 Pb. We have decomposed the contributions to the strength into OfkJ, 
and 2fvjj and higher components and found that the Ofiw contributes mainly to the resonance. 
The large number of 2hu and higher transitions form a very broad structure at excitation 
energies between 30 and 60 MeV. The calculated energy centroids are in very good agreement 
with values previously obtained using two different Skyrme interactions. The dependence of 
the centroids on neutron number has been shown for the Sn isotopic chain with the neutron- 
to-proton ratio in the interval from N/Z = 1 to N/Z = 2. The IVSM centroids are located 
at high excitation energies, ranging from 25 MeV above the GT centroids for the lightest 
isotopes, to 12 MeV for the heaviest. 

Evaluations of GT strength from experimental cross sections of charge-exchange reactions 
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take into account the isovector spin monopole mode but, because of the unknown distribution 
of the IVSM strength, its contribution is subtracted incoherently from the total measured 
strength. To analyze the validity of this procedure, we have calculated the L = strength 
using the GT+IVSM operator Eq. (fl6l) and the full L = operator Eq. (fl4"j) . It has been 
found that the inclusion of the isovector spin monopole term contributes to the strength 
at high excitation energies, and also reduces the strength of the resonance. The shift of 
the strength to higher excitation energies has been analyzed for the Sn isotopic chain with 
masses in the range 100 < A < 150. The total L = strength for isotopes with low number 
of excess neutrons is enhanced, whereas it is not modified for isotopes with A > 120. The 
full L = operator only changes the strength at high excitation energies, i.e. for large 
momentum transfer. A similar analysis of the effect of finite momentum transfer has been 
performed for 48 Ca, 90 Zr and 208 Pb. For 90 Zr, in particular, the results have been compared 
with a recent analysis of the GT quenching factor based on (p, n) and (n, p) data. We have 
found that the total strength in the /3~ channel is not modified by the inclusion of higher 
order terms, i.e. S9J = S^T°. The strength in the /3+ channel is reduced by approximately 
15% by the inclusion of finite momentum-transfer terms, contrary to the assumption made 
in the analysis of experimental cross sections. Combining these results with the model 
prediction that 6% of the strength in the j3~ is located above 50 MeV excitation energy, we 
find that the Ikeda sum rule is satisfied within experimental uncertainty. 

The determination of Gamow- Teller strength, complicated however by the excitation of 



the IVSM mode, was also performed using the ( 3 He,t) reaction on 208 Pb [34] and 150 Nd 35]. 
Because the IVSM transition density has a node close to the surface, probes that penetrate 
deep into the nucleus display smaller cross sections due to the cancellation of contributions 
from the surface and the bulk. In contrast, probes absorbed at the surface have larger cross 
sections because there is no contribution from the volume region. An analogous effect can 



be obtained using probes with different energies (see Sec. IV. A in Ref. [1J|). The present 
calculation does not differentiate between various probes, and the only effect of the energy 
of the incoming probe is on momentum transfer. It would be interesting to perform a study 
of the interference of Gamow- Teller and isovector spin-monopole modes, taking into account 
the characteristics of the experimental probe. One could, in particular, combine the RQRPA 
transition densities with a DWIA calculation, and compare the resulting cross sections with 
the experiment. 
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